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Abstract. I will present a review of the field of Hanbury Brown- Twiss 
interferometry in relativistic heavy-ion collisions. The "HBT puzzle" is explored 
in detail, emphasizing recent theoretical attempts to understand the persisting 
puzzle. I also present recent experimental results on azimuthally sensitive HBT, 
HBT of direct photons, and some surprises in the comparison of HBT results from 
p+p and Au-f Au collisions at RHIC. 



PACS numbers: 25.75.Gz 
1. Introduction 

Among the large set of measurements predicted ^ to uncover the formation and 
nature of a quark-gluon plasma (QGP) in relativistic heavy-ion collisions, and 
the corresponding phase transition between this plasma and the more familiar 
hadronic phase, Hanbury Brown-Twiss interferometry (HBT) |2] surely is the most 
maligned. This Quark Matter conference marks the three-year anniversary of the 
first presentation of experimental HBT results from RHIC, which showed that 
the "standard" HBT source radii in Au-fAu collisions at ^/s = 130 GeV (introduced 
below) arc quantitatively similar to previous measurements at both AGS and SPS. 
The CERES Collaboration illustrated this nicely in results 0] presented at the last 
Quark Matter. The lack of an energy dependence, the subsequent inconsistency 
with otherwise-successful hydrodynamical models, and the failure of most alternative 
explanations, led to the moniker "HBT puzzle" 

Whether the "HBT puzzle" continues to be cause for concern in our field remains 
to be seen, as it has recently taken a back seat to studies that do show strong differences 
at RHIC compared to lower energies (most notably higher px effects attributed to 
jet quenching [7]). The puzzle has also been somewhat cast aside as a simple lack 
of adequate description for the freeze-out process (most notably at this conference 
by Miklos Gyulassy, an early advocate of HBT as phase transition signature |H1). 
While it is clear that understanding the evolution of the hot, dense, expanding system 
requires an understanding of its space-time structure accessed by HBT studies, it is not 
clear to me whether this understanding can be bypassed when making QGP discovery 
claims based on other obscrvablcs. In that sense, the HBT puzzle still warrants much 
theoretical and experimental effort, perhaps more than it currently receives. 

Providing a ray of optimism, Rcinhard Stock, in his Quark Matter introductory 
overview, characterized HBT as bringing us "from enlightenment to tragedy, and back 
to enlightenment, on a time scale of four years." With these words in mind, in this 
overview I assess the status of HBT measurements and the persisting puzzle, mostly 
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Table 1. Some of the HBT differential studies underway or undertaken recently 
in hoavy-ion collisions. 



Diff. quantity 


What it investigates 


Recent results 


Beam energy 


Onset effects, transition phenomena 


Him 


Tranverse momentum 


Dynamics, collective expansion 




Particle type 


Hydrodynamic my scaling 




Collision system 


Origin of Bose-Einstein enhancement 


rrei 


Azimuthal angle 


Spatial anisotropy, system evolution 


EH 



from an experimentalist's standpoint. I also discuss new developments in two-particle 
interferometry that contribute more to our understanding of heavy-ion collisions at 
relativistic energies. (For an experimental summary of related results presented at 
this conference, see H. Appelshauser's contribution to these proceedings [S].) 

2. HBT in relativistic heavy-ion collisions 

The experimental technique of using two-particle interferometry to relate the 
momentum space separation of particles to their separation in space-time is well- 
established In the case of identical bosons, e.g. tt"*" mesons, quantum interference 
among the particles leads to an enhancement of pairs with small momentum difference 
q (Bose-Einstein enhancement). To isolate the small set of correlated pairs that 
undergo this quantum interference from the enormous amount of uncorrelated pairs 
in an event, a correlation function C(q) is formed in which pairs from real events 
are divided by pairs from different events. In heavy-ion collisions, C(q) is often 
constructed in three dimensions and fit to a three-dimensional Gaussian: 

^, , real pairs ,^r, , „2 p2 2 „2 2 p2 , 

C(q) = — f- — — = 1 4- Ae"«°"'^^°"'"'?-dc«=idc-9io„g«io„g ^ (1) 

mixed pairs 

where the subscripts indicate the long (parallel to beam), side (perpendicular to beam 
and total pair momentum k) and out (perpendicular to qi and qs) decomposition of q. 
A'^ is a normalization constant. The i?'s in Eq. ^ known as the HBT radii, quantify 
the widths of the Gaussians and represent the apparent size of the particle source, 
which may depend on the transverse momentum slice under study {i.e. homogeneity 
regions UOI). In practice, final-state effects such as Coulomb also contribute to C(q) 
and need to be accounted for. I won't discuss these here except to note that nearly 
all heavy-ion studies now have adopted an improved Coulomb treatment 

The purpose of HBT studies in heavy-ion collisions is to explore the space- 
time evolution and freeze-out of the system. This can be thought of as three-fold: 
the spatial distribution of the emission points, the time length of emission, and 
the dynamical properties of the system as it evolves. HBT serves as a tool for 
disentangling these contributions, and the out-side-long decomposition of q is chosen 
for that reason. Experimentally, HBT radii are studied as differentially as statistics 
and detector configurations allow; see Table 1. At RHIC, greatly increased pion 
production [dN/dy ~ 300 per flavor ^H]) due to higher collision energies makes triple- 
differential HBT analyses possible, e.g. kx x ccntrality x (f>paii- |17| . 

Highlighting two differential studies in particular: (a) The transverse momentum 
(fcr) dependence of the HBT radii for identical pions probably is studied most 
often, under the model-dependent view that space-momentum correlations in the 
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Figure 1. (a) Contours representing densities of omissions points in tlie 
transverse plane for kx = 0.0 GeV/c (top) and fcy = 0.5 GcV/c (bottom). Tlie 
emission direction is to the right, illustrating that higher px particles emerge 
nearer to the surface, (b) Contours of constant energy density at two different 
times in the evolution of a noncentral collision. Both figures are hydrodynamic 
calculations taken from Ref. 1191 . 

source are due mostly to collective expansion |2|- As the source expands, radial 
flow pushes higher px particles more at the surface (Figure 1(a)). Within this 
picture, analytical expressions have been derived to extract the expansion velocity 
and emission duration from the mx (jut = \/p'^ + rn?) dependence of the HBT 
radii. This is discussed further in section 3. (b) HBT studies relative to the reaction 
plane in non-central collisions allow the possibility to compare the expanded system's 
transverse eccentricity at freeze-out to its initial eccentricity from a nuclear overlap 
model calculation (Figure 1(b)). This is discussed further in section 4. 

3. The HBT Puzzle 

The hydrodynamical approach to understanding HBT is motivated at RHIC by the 
model's demonstrated ability to describe soft pT spectra and elliptic flow consistently 
for several particle species ^Hl- "Hydro" calculations for these observables point to 
fast thermalization in a partonic phase, followed by hydrodynamic expansion for ^ 15 
fm/c with an intermediate phase transition. However, these calculations yield strong 
disagreement with HBT radii i?out and i?iong are overpredictcd by as much as 
a factor 2, and i?side is somewhat underpredicted. In particular, the measured fc^ 
dependence of -Rside is in contrast to hydro and other models that predict little (if 
any) fc-r dependence. This disagreement, and the lack of energy dependence of the 
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HBT radii for a fixed bin is known as the "HBT Puzzle." Here I will summarize 
recent theoretical approaches to resolving the puzzle, illustrating along the way why 
the data are indeed puzzling. 

The first question is whether the hydro calculation itself can be altered to agree 
better with the HBT radii. This was explored in great detail by Heinz and Kolb |19j . 
A summary of their findings: 

• Default initial conditions^ such that pt spectra and elliptic flow arc well-described, 
yield the disagreement discussed above. 

• Freeze- out directly at hadronization brings i?out and i?iong close to the data, but 
at the cost of large disagreement with pr spectra. -Rsidc doesn't move, still lacking 
fcr dependence. 

• Faster thermalization or non-zero initial flow also reduces i?out ^md i?iong 
compared to default conditions, but not enough to agree with measurement. 

Hirano and Tsuda \2,\\ also checked the effect of maintaining proper particle 
abundances with chemical potentials in the hadronic phase of hydrodynamic evolution. 
They found better agreement for i?iong and i?out, but i?side disagreed more. 

Whether or not hydro can be made to reproduce HBT radii while maintaining 
the strong agreement with momentum-space quantities might come down to the 
assumption of longitudinal boost invariance implicit in most hydro approaches. Csorgo 
has shown [221 that, by introducing a Hubblc-like flow and allowing for a smeared 
freeze-out temperature, the HBT radii (as well as spectra and elliptic flow) can be 
well described by a single parameter set. However, a blast-wave parametrization pS) 
which maintains longitudinal boost invariance also was able to fit the data rather 
well. Both of these "hydro-inspired" approaches do not contain a full hydrodynamic 
evolution, but the descriptions of freeze-out obtained by their fits possibly hint at the 
directions full hydro models should take. 

A few alternatives to hydro have been investigated. Studies of the effect of 
introducing opacity in a parton cascade model |24| showed that the pion freeze-out 
distribution is indeed sensitive to the transport opacity in the partonic phase. For a 
parton cascade with no opacity, all three HBT radii show no dependence and values 
below the data. As the opacity is increased, -Riong and i?out develop fcy dependences 
and increase toward the data nicely, but i?sidc doesn't move. (The authors of Ref. [21] 
claim this points to a lack of sensitivity of i?sido to early partonic dynamics.) Still, 
the behavior of i?iong and i?out with increasing opacity may indicate that reality lies 
somewhere between the extremes of cascade and ideal hydro. 

Finally, a natural consequence of hydro models is a negative x-t correlation, i. e. 
pions further from the source's center are emitted earlier. i?out depends explicitly on 
an [xt) cross-term, with a negative {xt) acting to increase i?out- However, cascade 
models such as AMPT |21] have shown that positive x-t correlations arise naturally 
in their codes, thereby reducing i?out compared to hydro's overprediction. Whether 
it makes physical sense that particles emitted further from the center could decouple 
after particles emitted closer to the center is an important question that needs to be 
confronted by these models. 

4. Recent advances 

Despite our difficulties to understand what standard pion HBT measurements are 
telling us at RHIC (and therefore at lower energies, for that matter), two-particle 
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Figure 2. Source eccentricity obtained with azimuthally-scnsitivc HBT (cfinal) 
vs. initial eccentricity from a Glauber model (einitial) in Au+Au collisions at y/s = 
200 GeV. The most peripheral collisions correspond to the largest eccentricity. 
The dashed line indicates Cinitiai — ^finai- 

IT7I . 



intcrfcromctry studies have been extended recently into new domains at both SPS 
and RHIC. Here I wiU discuss two examples. 

STAR has recently completed an analysis of the azimuthal dependence of HBT 
radii relative to the reaction plane at y/s^N = 200 GeV ^7|. Azimuthally sensitive 
HBT was suggested as a probe of how spatial anisotropy evolves in non-central 
collisions. The reasoning is straightforward: (a) The initial almond-shaped geometry 
gives rise to anisotropics in pressure gradients, the same gradients responsible for 
elliptic flow, (b) The pressure gradients drive a preferential expansion in the reaction 
plane that decreases the spatial anisotropy. (c) HBT provides a measure of the freeze- 
out source shape, which in principle could change its orientation from out-of-plane to 
in-plane extended depending on the amount of pressure built up and the expansion 
time. 

star's results showed an intuitive (though model-dependent) centrality 
dependence of the system's eccentricity at freeze-out (Figure|2J). Near-central collisions 
showed final eccentricities consistent with zero. When going from central to peripheral 
collisions, the final eccentricity increased, reflecting the greater initial eccentricity 
while retaining some of its initial out-of-plane almond shape. Given the strong evidence 
for significant pressure build-up in the system from elliptic flow measurements, the 
results point to short evolution times as the dominant cause for out-of-plane freeze-out 
shapes. 

Another domain that holds much promise is HBT of direct photons Unlike 
pions, 7-HBT probes the initial state, thereby providing a potential probe of conditions 
in the deconflned phase. However it is a major challenge experimentally, due to the 
low relative dircct-7 yield and the many sources of pairs with small relative momenta. 
To name a few: photon conversions, decays of tt'^'s (which themselves are affected by 
HBT), resonance decays, misidentified 7's. In spite of the technical challenges, WA98 
recently performed \28\ the first 7-HBT measurement in heavy-ion collisions, in a low- 
momentum range of 0.15 < pr < 0.35 GeV/c. Very interestingly, they observed HBT 
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radii quantitatively similar to tt-HBT for the same pt, which they attributed to soft 
photons arising in late stages of the collision. 

Though this was surely already significant enough for a Phys. Rev. Lett., WA98 
took the analysis a step further, using the A parameter from fitting the correlation 
function (and assuming a fully chaotic source) to determine the absolute direct photon 
yield in this momentum range. This technique provides a complementary tool to 
subtraction-based direct-7 yield measurements, and at higher transverse momenta 
it holds great promise for accessing the temperature reached at early stages of the 
collision. (Direct photon studies are underway at RHIC.) 

5. HBT and freeze-out from p+p to Au+Au 

One of the advantages of the RHIC experimental program is the ability to collide 
different systems with the same center-of-mass energy, allowing for identical analyses 
of these different systems. Preliminary results were presented at this conference jl5) 
of pion HBT in p+p and d+Au collisions at y^ijvw = 200 GeV, for comparison to 
the centrality-binncd Au+Au analysis (with fixed particle id. cuts, same kr bins). 
Like the heavy-ion case, the three HBT radii in p+p exhibit a characteristic decrease 
with increasing kx, though for p+p it has been attributed to string and multistring 
fragmentation in earlier studies. A somewhat surprising result comes about when 
dividing the Au+Au and d+Au radii by the p+p radii ^Sl- The divided trends are 
roughly flat with kx for all radii, indicating an apparent scaling in the /ct dependences 
of the HBT radii for these three systems. Given that the kx dependence presumably 
arises in very different ways (collective expansion in heavy ions, string fragmentation 
in elementary particles), these results also are a bit puzzling. 

Similarity among these systems can be further tested. At the last Quark Matter 
conference, CERES proposed a simple ansatz to investigate if pion freeze-out is 
driven by a critical mean free path A/. They developed an approximate expression 
for A/ that depends on the freeze-out volume Vf estimated from HBT radii, pion and 
nucleon particle densities (Ni), and scattering cross sections for these particles (cry ). 
Here I only give the final expression for A/; see their paper j^H] for full development 
of, and reasoning behind, this expression and its caveats (such as how the Ni are 
estimated from the dN/dyi, i?iong5 a-nd thermal velocity /3th): 

^ ^Vf ^ (2^)^/^i?longi^Le (2) 

Figure ini^a) shows the numerator and denominator of Eq. [3 plotted in the same panel 
vs. ^/s (from AGS to RHIC energies), with the scale for Vf {Na) defined on the 
right (left) axis. Both trends exhibit non-monotonic behaviors with dips in between 
AGS and SPS energies that match impressively for the given scales. Vf decreases at 
AGS energies mostly due to decreasing i?sidc and increases from SPS to RHIC due to 
increasing i?iong- The Na expression first decreases and then increases, reflecting the 
transition between nucleon and pion dominance as the chemical composition changes 
with increasing yG. The agreement between the trends yields a critical mean free 
path of ~ 1 fm across the energy range, which the authors attribute to significant 
opaqueness in the source [221 ■ 

Using the HBT radii for p+p and d+Au presented at this conference and 
published dN/dy values from STAR [30], the critical mean free path ansatz can be 
applied to these systems. Figure |2b) shows the incredible result that A/ is ~ 1 
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Figure 3. (a) The freeze-out volume Vf and the density /cross-section term Na 
as a function of y/s for central heavy-ion collisions, (b) The same quantities 
calculated for p+p and d+Au at y/s = 200 GeV, plotted on the same relative 
scale; circles (•) are Na, stars (*) are Vf. 



fm for these light systems as well, even though Vf and Na arc lower by two orders 
of magnitude. I find this really incredible. Why would these systems exhibit such 
similar mean free paths at freeze-out in this simple ansatz? Is it only a coincidence, 
for example that A/ is reflecting the considerable opaqueness in Au-|-Au but merely 
the system size in p+p and d-f Au? Or has CERES touched upon a deeper connection 
between the freeze-out volume and freeze-out density that is somehow independent of 
the actual dynamics in the system? 



6. Final remarks 



The status of HBT interfcrometry as a signature for the formation of a quark-gluon 
plasma or the corresponding phase transition remains murky. This murkincss is due 
mostly to the persisting disagreement between tt-HBT data and most hydrodynamical 
models that obtain good agreement with transverse momentum spectra and elliptic 
flow, as well as HBT's intrinsic sensitivity to the latest stages of the collision. After 
discussing the HBT puzzle ad infinitum with many experts before and during this 
conference, I am of the opinion that there are primarily two candidate solutions: (a) 
the growing evidence for a lack of longitudinal boost invariance at RHIC suggests that 
alternative hydro formulations need to be explored further; (b) parton cascade studies 
that are able to get close to the HBT data with finite opacity might indicate that, 
contrary to popular belief, perhaps the hydrodynamic "limit" has not been reached at 
RHIC. Both statements illustrate the importance of HBT in understanding heavy-ion 
collisions. 

However, a number of recent developments have shed light on HBT's ability to 
unlock information about the nature of the collisions we study. Aziniuthally-sensitive 
HBT measurements of non-central collisions indicate that the system retains some 
of its initial almond shape at the end of its evolution, supporting the notion that 
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the system lives shorter than predicted. First observation of direct photon HBT 
at SPS energies illustrates the potential for 7-HBT to probe early conditions in the 
collision. Identical experimental conditions and analysis methods for p+p and Au+Au 
at RHIC provide an important link between these two systems with presumably very 
different dynamics. The larger Au+Au dataset taken this year at RHIC also holds 
the potential for exotic correlation studies between non-identical particles, further 
and more differential measurements of "standard" HBT extended to higher transverse 
momenta, and maybe even studies of the quark gluon plasma using HBT of direct 
photons. Though HBT interferometry has been troublesome to many of us in the first 
years of RHIC running, its future contributions to our understanding of relativistic 
heavy-ion collisions are assured. 

Many people helped in the preparation of this overview. In particular I would like 
to thank S. Bekele, T. Gutierrez, M. Heffner, B. Holzmann, and A. Kisiel for supplying 
figures on short notice. I especially want to thank H. Appclshauscr, T. Csorgo, 
U. Heinz, A. Kisiel, M. Lisa, D. Molnar, S. Soff, and D. Teancy for enlightening 
discussions. This work was supported by NSF grant PHY-0099476. 
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